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Abstract 

Nuclear fusion is a promising energy source of the future. One possible 
way to achieve reliable energy production from fusion reactions is to 
confine the high temperature deuterium-tritium fuel by magnetic fields. In 
magnetic confinement fusion devices like tokamaks, it is crucial to confine 
the high energy fusion-born helium nuclei (a-particles) to maintain the 
energy equilibrium of the plasma. However, super-thermal energetic ions can 
excite various instabilities which can lead to their enhanced radial transport. 
Consequently, these instabilities may degrade the heating efficiency and they 
can also cause harmful power loads on the plasma-facing components of the 
device. Therefore, the understanding of these modes is a key issue regarding 
future burning plasma experiments. 

One of the main open questions concerning energetic particle (EP) driven 
instabilities is the non-linear evolution of the mode structure. In this thesis, 
I present my results on the investigation of /^-induced Alfven eigenmodes 
(BAEs) and EP-driven geodesic acoustic modes (EGAMs) observed in 
the ramp-up phase of off-axis neutral beam injection heated plasmas in 
the ASDEX Upgrade tokamak. These modes were well visible on several 
line-of-sights of the soft X-ray cameras which made it possible to analyse the 
time evolution of their spatial structure. 

In order to investigate the radial structure, the mode amplitude has to 
be determined on different line-of-sights. I developed an advanced ampli¬ 
tude reconstruction method which can handle the rapidly changing mode 
frequency and the low signal-to-noise ratio. This method is based on short 
time Fourier transform which is widely applied in the thesis, because it is 
ideal to investigate the time evolution of transient wave-like phenomena. 

The radial structure analysis showed that in case of the observed 
downward chirping BAEs the changes in the radial eigenfunction were smaller 
than the uncertainty of the measurement, while in case of rapidly upward 
chirping EGAMs the analysis shows shrinkage of the mode structure. These 
experimental results are shown to be consistent with the corresponding 
theory. 



Gyors ionok hajtotta hullamok vizsgalata tokamak plazmakban 

Horvath Laszlo 

Nuklearis Technikai Intezet 
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Kivonat 

Jovonk igeretes energiaforrasa a nuklearis fuzio, mely megvalositasanak 
egyik lehetseges modja magneses osszetartasu berendezes alkalmazasa. Ilyen 
berendezes a tokamak, melyben a deuterium-tricium fiizios reakcio soran 
keletkezo nagy energiaju helium atommagok (a-reszecskek) osszetartasa 
kiilonosen fontos a plazma energiaegyensulyanak fenntartasa szempontjabol. 
Ezek a termikusnal joval nagyobb energiaval rendelkezo ionok azonban olyan 
plazmainstabilitasokat kelthetnek, melyek felgyorsithatjak a gyors ionok 
radialis transzportjat. Ennek kovetkezmenyekent csokken az ct-reszecskek 
futesi hatekonysaga, tovabba a kiszorodo gyors ionok sulyosan karosfthatjak 
a plazmat koriilvevd falelcmeket. Emiatt a gyors ionok altal hajtott in- 
stabilitasok megertese letfontossagu a jovo energiatermelo fiizios reaktorai 
szempontjabol. 

A kutatasi tertilet egyik maig nyitott kerdese a gyors ionok altal 
hajtott modusok szerkezetenek nem-linearis idofejlodese. Jelen dolgozatban 
az ASDEX Upgrade tokamak semleges atomnyalab altal ftitott plazmainak 
kezdeti szakaszaban megfigyelt /Uhajtott Alfven sajatmodusokkal (/3-induced 
Alfven eigenmode, BAE) es gyors reszecskek altal hajtott geodezikus 
akusztikus modusokkal (energetic particle driven geodesic acoustic mode, 
EGAM) kapcsolatos kiserleti eredmenyeimct mutatom be. Ezen modusok 
jol megfigyelhetoek a lagy-rontgen diagnosztika szamos csatornajan, mely 
segitsegevcl a modusok radialis strukturajanak idofejlodese vizsgalhato. 

A radialis struktura tanulmanyozasahoz a modus amplitudo 
meghatarozasara van sztikseg a lagy-rontgen diagnosztika ktilonbozo 
latovonalain. Ennek erdekeben kifejlesztettem egy amplitudo rekonstrukcios 
cljarast, amcly segitsegevcl kezelhcto a modus gyorsan valtozo frekvenciaja es 
az alacsony jel-zaj viszony. Ez a modszer rovid idejii Fourier-transzformacion 
alapul, mclyet szeles korben hasznaltam a munkarn soran, mivcl kivalo 
tranziens jclck elcmzesere. 

Vizsgalataim megmutattak, hogy a megfigyelt BAE-k eseteben a modus 
radialis strukturajanak valtozasa kisebb a meres bizonytalansaganal. Ezzel 
szemben az EGAM-ok eseteben a modus struktura osszehuzodasa figyclheto 
meg. Ezek a kiserleti eredmenyek osszhangban vannak a jclcnseget leiro 
clmelettel. 



A szakdolgozat kifrasa 

Azonosfto: DM-2014-44 

Onfenntartoan mukodo szabalyozott magfuzios reaktor mukodese szem- 
pontjabol kulcsfontossagu a termikusnal joval nagyobb energiaju ionok 
osszetartasa. Ezen ionok rezonans reszecske-hullam kolcsonhatason keresztiil 
kiilonbozo hullamokat tudnak hajtani, mely hullamok a reszecskek meg- 
no vekedett veszteseget okozzak, ezaltal rontva a plazmafuzes hatekonysagat. 
A hallgato elsodleges feladata a nemetorszagi ASDEX Upgrade tokamakon 
megfigyelt, gyors reszecskekhez kotheto kiilonbozo Alfven-sajatmodusok, 
ezen beliil a gyors frekvenciavaltozassal jcllemezhetd un. „csdrpold” (az angol 
„chirp” szobol) modusok vizsgalata. Nyitott kerdes, hogy a csorpolo modusok 
radialis sajatfiiggvenye hogyan valtozik a csorp ideje alatt: a meresi bizonyta- 
lansagon beliil allando marad, vagy szignifikansan valtozik? Ez az informacio 
kulcsfontossagu a modusok megertese, es kesobbi prediktiv modcllczese szern- 
pontjabol. A hallgatonak ki kell dolgoznia egy eljarast amellyel a csorpolo 
modusok sajatfiiggvenyenek idobcli valtozasa vizsgalhato, es a valtozas / nem 
valtozas tenye a meresi bizonytalansagon beliil szignifikansan kimutathato. 
A vizsgalat soran elsodlegesen felhasznalando adatok az ASDEX Upgrade 
tokamakon rendelkezesre alio magneses diagnosztikak jclci, a vonalintegralt 
lagy rontgen sugarzas merese illetve az elektron-ciklotron kibocsajtason 
alapulo lokalizalt homersekletmeres. Sziikseg eseten mas diagnosztikak is 
bevonhatok a vizsgalatba. A fcladat soran meg kell ismerni a diagnosztikak 
miikodeset, ki kell dolgozni egy eljarast a radialis sajatfiiggveny idobcli 
valtozasanak meresere es a meresi bizonytalansagok becslesere. A meresi 
eredmenyek fiiggvenyeben, azok melyebb megertesehez sziiksegesse valhat 
a vizsgalt modusok szimulacioja es szintetikus diagnosztikak alkalmazasa. 
Ezen szimulacios programok megirasa a hallgatonak nem feladata, a celra jol 
miikodo programcsomagok allnak rendelkezesre a temavezeto intezeteben. 
Az elemzesek vegrehajtasahoz segitseget nyujt a BME NTI-ben fejlcsztett 
NTI Wavelet Tools programcsomag. A programcsomag tobb kiilonbozo 
clernzo algoritmus mellett jol miikodo kiserleti adatbeolvaso-, manipulalo- 
es abrazolo kepessegekkel rendelkezik, ezaltal celszerii a vizsgalathoz a mar 
meglevo programokat felhasznalni, sziikseg eseten azokat modosftani. Az 
NTI Wavelet Tools nagy resze IDL programnyelven irodott, ezert ennek a 
nyelvnek ismerete vagy sziikseges elofcltetel, vagy megtanulando, mivel a 
hallgato feladata lesz a sziikseges valtoztatasokat is elvegezni. A szakdolgozat 
temajahoz kotodo kutatomunka nemzetkozi egyiittmukodesben zajlik, eseti 
kiilfoldi kiutazasokkal. Magas szintii angol nyelvtudas, alapos programozoi 
es matematikai ismeretek sziiksegesek. 
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Chapter 1 
Introduction 


Our world provides plenty of options to meet the energy needs of mankind. 
One of them, which is already in the service of humanity, is the release of 
the binding energy of atomic nuclei. Energy is produced by nuclear fission 
in nuclear power plants all over the world. However, sustainable energy 
production to the electricity grid by fusion reactions has not yet been 
achieved [1]. 


1.1 Fusion power generation 


An example of reliable fusion energy production is the Sun. Its hydrogen ( 1 H) 
fuel is converted to helium ( 4 He) in different astrophysical reaction chains [2] . 
These fusion reactions provided energy for several billions of years to life 
on Earth. In stellar cores the rate of fusion reactions is very slow, because 
proton-neutron decays have to occur. Therefore, astrophysical reaction chains 
have very low cross-sections, thus the released power density of the Sun is 
less than 1 W/m 3 [3] . 

Reproduction of stellar core conditions on Earth for nuclear fusion power 
production would be inefficient. In order to achieve sufficiently high power 
from fusion reactions for electricity production, light isotopes should be 
chosen [3] such as deuterium ( 2 H = D), tritium ( 3 H = T) or 3 He. Some 
of the fusion reactions possible with the aforementioned isotopes are listed 
here [4]: 


D + D 3 He (0.82 MeV) + n (2.45 MeV), 

D + D => T (1.01 MeV) + p (3.02 MeV), 

D + T =>• 4 He (3.52 MeV) + n (14.1 MeV), 

D + 3 He =► 4 He (3.66 MeV) + p (14.6 MeV). 
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Considering the reaction rates presented in figure 1.1, the most favourable 
reaction on Earth is the deuterium-tritium (DT) reaction [5]. The reaction 
rate of DT peaks at a lower temperature and at a higher value than other 
reactions commonly considered for fusion energy. 



Figure 1.1: The fusion reaction rates of different isotopes [6]. 

A huge amount of deuterium is available in nature, because there is 1 atom 
of deuterium for every 6700 atoms of hydrogen in natural water resources [4], 
In addition to the enormous resources, deuterium can be easily extracted at 
a low cost. However, gathering tritium is difficult. Tritium is a radioactive 
isotope and it has a half-life of 12.3 years. Therefore, tritium can be found 
in nature only in trace quantities [7]. Given that there is no practical tritium 
breeding technology which can satisfy the demand of an energy producing 
fusion reactor, the required tritium must be produced by the reactor itself. 

Taking advantage of the fusion reaction generated neutrons, tritium can 
be produced from lithium using the following processes [4]: 

6 Li + n (thermal) = ^He + T 
7 Li + n (energetic) = 4 He + T + n 

The tritium breeding process is illustrated in figure 1.2. The fusion-born 
neutrons, by hitting the wall containing lithium, produce tritium. In this 
way the fuel of a fusion power plant is deuterium and lithium. Both are 
non-radioactive elements and are available in large quantities on Earth. The 
final product of the fusion reaction and the tritium breeding reaction are 
helium, thus neither the fuel, nor the end-product of the fusion and breeding 
cycle is radioactive. 





Figure 1.2: Operational scheme of tritium breeding. The fusion-born neutrons, 
by hitting the wall containing lithium, produce tritium. 


Fusion experiments with particle accelerators provided useful information 
to determine the cross sections of the aforementioned reactions. However, 
sustainable energy production with accelerators cannot be achieved [6], since 
the cross-section of Coulomb scattering is several orders of magnitude higher 
than the cross-section of the fusion reaction. Therefore, the most promising 
method is to heat the DT fuel to sufficiently high temperature that the 
thermal velocities are high enough to produce fusion reactions. In this case 
the Coulomb collisions just redistribute the energy between particles. The 
necessary temperature is in the order of hundred million Kelvin for DT fuel. 
The conditions required for a self-sustaining process are estimated by the 
Lawson criterion for DT fuel and optimum temperature (~25 keV) [6]: 

ute > 10 20 s/m 3 , (1.1) 

where n is the density and Te is the energy confinement time which measures 
the rate at which a system loses energy to its environment: 

W 

te = p-, (1-2) 

-Moss 

where W is the plasma energy content and Pi oss is the total power loss of the 
plasma. 
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1.2 Magnetic confinement fusion 

One way of confining the fuel which has such high temperature is to employ 
a magnetic field. At temperatures necessary for the fusion reaction, the 
particles are ionized and the fuel is in the plasma state. The charged particles 
(ions and electrons) start their helical motion in the plane along the magnetic 
field. The circular part of this motion in the plane perpendicular to the field 
line is called gyromotion. The helical motion is suitable to confine the plasma, 
because it does not allow the particles to move perpendicular to the magnetic 
field, but they are free to move parallel to the magnetic field. The problem of 
losses at the ends of a linear device can be solved by bending the magnetic 
field lines into a torus. However, in a toroidal geometry, the magnetic field is 
not homogeneous. The strength of the magnetic field is inversely proportional 
to the distance measured from the centre of the torus, because of the higher 
current density in the toroidal held coils (indicated with red in figure 1.3) 
in the inner side of the torus results in a stronger magnetic held there. This 
inhomogeneity leads to a charge dependent drift (VH drift) [6], which causes 
the electrons and ions to move in vertically opposite directions. The resulting 
vertical electric held creates a charge independent E x B drift [6] , which moves 
the entire plasma towards the outside of the torus. This electric held between 
the top and the bottom of the torus can be shorted out by helically twisting 
the magnetic held lines. 

A type of magnetic confinement device is the tokamak [6], where the 
magnetic held lines are helically twisted by a plasma current as it is illustrated 
in hgure 1.3. In tokamak geometry the toroidal direction is the long way, the 
poloidal direction is the short way around the torus. The toroidal plasma 
current is driven by a transformer coil and it twists the magnetic held 
lines because it generates a poloidal magnetic held indicated with purple in 
hgure 1.3. The twist of magnetic held lines is described by the safety factor. 
The safety factor q is the ratio of toroidal transits per single poloidal transit 
of a magnetic held line [8]. Due to the axisymmetry of a tokamak device, 
the magnetic held lines are organized into magnetic hux surfaces. Along the 
magnetic held lines, the particles can move freely, thus the plasma parameters 
on the magnetic hux surfaces are balanced very quickly, while the transport 
perpendicular to the surfaces is several orders of magnitude slower. 

In order to realize a self-sustaining magnetically confined plasma, the 
released fusion energy has to heat the plasma. The energy generated in 
a fusion reaction is distributed as kinetic energy to the fusion products. 
According to their mass ratio, 80 % of the energy is carried by the neutrons 
which do not interact with the plasma and hit the plasma facing components. 
The remaining 20 % energy is carried by the He nuclei. These He 2+ ions have 
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Figure 1.3: Schematic drawing of a tokamak. The superposition of the toroidal 
magnetic field and the poloidal magnetic field generated by the 
toroidal plasma current leads to helically twisted magnetic field lines 
in the torus. 


3.5 MeV initial energy and isotropic angular distribution. The self-heating 
provided by the He 2+ ions (a-particles) is called a-heating. In a hot plasma 
the typical mean free path of fusion-born a-particles is very long compared 
to the size of the device [9]. Thus, it is crucial to well confine the high energy 
a-particles. 

Beside the a-heating, auxiliary heating systems are needed to reach fusion 
relevant conditions. Furthermore, these external heating methods are planned 
to be used to tailor the plasma profiles during the fusion burn. Two types of 
auxiliary heating systems are used to generate fast ions, namely the neutral 
beam injection (NBI) and the ion cyclotron resonance heating (ICRH). These 
techniques are capable of accelerating hydrogen isotopes up to the MeV 
energy range. Since plasma heating scenarios are predominantly based on 
fast ions, their transport properties are fundamental to the success of any 
burning plasma experiment. Enhanced radial transport of energetic particles 
(EPs) degrades the heating efficiency and it can also cause harmful power 
loads upon the Erst wall of the device. Since, in the plasma core, the most 
important transport process of EPs is their interaction with global plasma 
waves [9], the thorough understanding of the EP driven instabilities is crucial. 
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1.3 Fast ion transport due to plasma waves 


Marginally stable plasma modes can be destabilized by EPs in tokamaks. In 
addition to normal modes, there are also energetic particle modes (EPMs) 
which are forced oscillations characterized by strong dependence on the fast 
ion distribution function. First, a basic introduction of the shear Alfven waves 
is presented. 

Shear Alfven waves are electromagnetic waves that propagate along the 
magnetic field with the following dispersion relation [10]: 

a; 2 = k\v\, (1.3) 

where 

Va = / y£ ( L4 ) 

V ^0 / j i TliTTli 

is the Alfven velocity, ku is the wave number in the direction of the magnetic 
field, n i m i i s the mass density of the plasma, i denotes the different ion 
species, B is the magnitude of the magnetic field and /iq is the vacuum 
permeability. In a torus, the periodic boundary condition requires that the 
wave solutions are quantized in toroidal and poloidal directions: 


k\ 


n — m/q{r ) 

Rn 


(1.5) 


where n and m are the toroidal and poloidal mode numbers respectively, q is 
the safety factor and i? 0 is the major radius of the plasma. In a toroidal 
plasma both k\\ and va are functions of the minor radius r. This leads 
to the so called shear Alfven continuum. The dispersion relation of two 
successive modes with poloidal mode number rn — 2 and 3 is illustrated on 
figure 1.4. Waves at different radii have different phase velocities which leads 
to phase-mixing and a hypothetical wave packet would rapidly disperse. This 
process is called continuum damping [11], In a torus with a large aspect ratio 
(i?o ^ a) the magnetic field can be approximated by £> ~ B 0 (l—e cos(6>)) [6], 
where e = a/Ro is the inverse aspect ratio and 6 is the poloidal coordinate. 
The parallel component of the wave vector now also depends on the poloidal 
coordinate 6 which can lead to intersection of the dispersion relations 
belonging to different poloidal mode numbers m. The crossing is expected 
to be at the radius where the wave number of the two counter-propagating 
waves equals: 

^mi(A) (1-6) 

The crossing of m — 2 curve (blue) and the rn — 3 curve (green) is visible 
in figure 1.4. However, toroidicity resolves this degeneracy at the points 
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Alfven continuum 



Minor radius [m] 

Figure 1.4: Shear Alfven continuum plotted (without coupling) with dashed blue 
line (n = 1, m = 2) and with dotted green line (n = 2, m = 3). 
In toroidal geometry, a frequency gap is formed due to the toroidal 
coupling between m = 2 and m = 3 modes. The Alfven continuum 
with coupling is plotted with solid red lines. In order to calculate these 
dispersion relations, the background plasma parameters are adopted 
from Lesur’s PhD thesis [11]. The major radius is 3.3 m, the minor 
radius is 0.96 m and the magnetic field is 1.2 T on axis. The q profile 
has been modeled by q(r) = 1.2 + 2.1(r/a) 1,5 and the electron density 
by n e (r ) = 0.11 + 1.57(1 — r 2 /a 2 ) 0 ' 3 [10 19 m~ 3 ]. 

of intersection and produces gaps in the frequency continuum as shown 
in figure 1.4 with red. For a particular toroidal mode number, and to the 
lowest order in the inverse aspect ratio (e = a/Ro ) expansion parameter, 
each poloidal harmonic couples to its nearest neighbouring sidebands [12] 
and the eq. (1.3) dispersion relations changes: 

_ _ k \m + k \\m +1 ± J( k \\m ~ k \\m+l) 2 + ^ k \\m k \\m+l 

“'A =-2(137)-■ ( Lr > 

Within these gaps global modes can exist, because waves do not experience 
heavy continuum damping. These modes are called toroidicity-induced Alfven 
eigenmodes (TAEs). Analogous to the TAE gap further gaps arise caused by 
the coupling of the m, m + 2 and m, m + 3 modes such as clipticity induced 
gap and non-up-down-symmetric gap [9]. In these gaps radially extended, 
weakly damped modes can exist, which can be driven unstable by the fast 
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ion population. The properties of these eigenmodes such as frequency or mode 
structure are mostly determined by the background plasma. A similar gap 
mode is the ^-induced Alfven eigenmode (BAE). The BAE gap is introduced 
by the coupling between compressible acoustic waves and the shear Alfven 
continuum [13]. 

These modes can be driven unstable by the free energy in the fast ion 
distribution function [9]. Energy transfer between an ion and a wave requires 
that the velocity vector is not perpendicular to the electric field of the wave. 
Thus, only the drift velocity of ions contributes to the energy transfer [14]. 
The energy transfer associated with the gyromotion averages to zero, since 
the gyromotion is very rapid compared to the mode frequency. 

The growth rate of a mode is generally small compared to the wave 
frequency, because significant net energy transfer is gained only after dozens 
of orbital cycles. To avoid phase averaging to zero, a harmonic of the 
drift-orbit frequency must match the wave frequency, so the following 
resonance condition has to be fulhlled for net energy gain [14]: 

u + (m + l)ue — nuj£ ~ 0, (1.8) 

where u> is the wave frequency, uJ q is the frequency of the poloidal motion of 
the ion, cis the frequency of the toroidal motion of the ion, m is the poloidal 
mode number and n is the toroidal mode number of the mode. I denotes 
the order of the harmonic of the drift-orbit frequency. In strongly shaped 
plasmas and for large drift-orbit displacements, higher order harmonics are 
also important. 

The net energy exchange between the wave and the particles is determined 
by several factors. It depends on the amplitude of the particular poloidal 
harmonic and on the drift velocity. Furthermore, the magnitude of the 
energy exchange also depends on the gradient associated with the particular 
resonance in the EP distribution function. Particles with speed slightly below 
the resonant velocity gain energy, while particles with speed slightly above the 
resonant velocity lose energy. Thus, EPs can drive a mode where the gradient 
in the EP distribution is positive along the coordinate defining the resonance. 
Mostly the gradient in the EP distribution respect to the energy is negative, 
so the energy distribution usually damps the wave. This is true even for the 
slowing-down distribution in a burning plasma. Modes are usually driven by 
the free energy in the spatial gradient. However, the anisotropy in the velocity 
space can also drive the mode. Finally, the alignment of the particle orbit 
and the eigenmode is also a governing factor of the energy exchange. Modes 
are excited if the drive can overcome the damping of the background plasma. 
Important damping mechanisms include ion Landau damping, continuum 
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damping and radiative damping [9]. The damping rates are very sensitive to 
plasma parameters such as the temperature and the q profile. 

Several types of non-linear behaviour of the mode amplitude and fre¬ 
quency are observed on present-day tokamaks. This behaviour ranges from 
steady-state - where the amplitude saturates at nearly constant frequency 
- to an explosive growth. The type of behaviour is expected to significantly 
influence the impact of the instabilities on the fast particle transport. Thus, 
their thorough understanding from both experimental and theoretical side is 
essential. In this thesis modes with bursting amplitude and rapidly changing 
mode frequency are experimentally investigated in detail. These modes are 
also called “chirping” modes. 

So far, the discussion was restricted to Alfven modes, but if the EP 
pressure is sufficiently high compared to the thermal pressure, the EP 
drive can overcome the continuum damping and energetic particle modes 
are excited. Their mode structure is independent from the shear Alfven 
continuum, and sensitively depends on the EP distribution. Their frequency 
usually correspond to a characteristic frequency of the EP orbital motion. 

One of the main open questions concerning EP driven Alfven eigenmodes 
and EPMs is the non-linear evolution of the mode amplitude. These instabil¬ 
ities constitute a non-linear system where kinetic and MHD non-linearities 
can both be important making it difficult to describe the phenomenon. In 
this thesis I present my results on the experimental investigation of beta 
induced Alfven eigenmodes (BAEs) and EP-driven geodesic acoustic modes 
(EGAMs) observed in the ramp-up phase of off-axis NBI heated plasmas in 
ASDEX Upgrade (AUG) [15]. 

The plasma scenario in which these modes were observed needs low 
density operation and off-axis NBI heating. Low density is important in 
order to keep the ratio of fast ion pressure and thermal pressure as high 
as possible. The excitation of these modes is possible if the damping is low 
enough. The damping of EGAMs is lower at lower temperatures. Thus, they 
are easily driven unstable at lower temperatures e.g. during the ramp-up 
phase of a discharge. The main damping mechanism of BAE is the ion Landau 
damping [16] which scales with the background ion temperature and density. 
This way in low density, low temperature plasmas BAEs are more easily 
excited. On AUG, even in low density operation, the damping of BAEs in 
the flat-top - where the main plasma parameters are maintained at constant 
values - is not sufficiently small to excite them with NBI. However, in the 
ramp-up phase BAEs are routinely driven unstable with off-axis NBI heating. 
The role of off-axis NBI is twofold. One is that, further out from the magnetic 
axis, the temperature and density is lower which again allows the excitation 
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of these modes. Second, the dominant auxiliary heating system in future 
burning plasma experiments will be off-axis NBI. 

The confinement of fast ions is crucial in a burning plasma experiment 
such as ITER where the heating of the plasma is based on the dominant 
self-heating by fusion-born o-particles. ITER is a large-scale scientific 
experiment under construction that aims to demonstrate the technological 
and scientific feasibility of fusion energy [17]. In general, EP driven modes 
have a massive contribution to the radial transport of fast ions. However, in 
view of that bursting BAEs and EGAMs are only observed in the ramp-up 
phase of discharges where damping is low enough due to the low density and 
temperature the question arises what are their implications regarding ITER 
or any other burning plasma experiment? 

BAEs are driven by the radial gradient in the EP distribution func¬ 
tion [16]. Thus, their presence enhances the radial transport of fast ions. 
However, their stability in ITER is not yet clear. Simulations have shown 
that in the ITER inductive baseline scenario the radial gradients in the 
EP distribution will not be sufficiently high to destabilize BAEs, but the 
difference is marginal and the uncertainty of these results is high [18]. EGAMs 
themselves do not cause radial redistribution of the fast ions directly, because 
they are driven by the velocity phase space gradient in the EP distribution. 
Therefore, the isotropic distribution of fusion-born a particles cannot affect 
EGAMs. However, NBI driven EGAMs can influence fast ion losses indirectly 
via mode-mode coupling with TAEs [19]. The impact of BAEs and EGAMs on 
the fast ion conhnement in ITER is not completely clear. On the other hand, 
the investigation of these modes is essential to understand the non-linear 
behaviour of EP driven modes. 

The main goal of this thesis is to experimentally investigate the rapid 
changes in the radial structure of bursting EP-driven modes during the 
non-linear chirping phase. Due to diagnostic and data analysis complexities 
this task has never been accomplished before. Some modes are expected to 
retain their radial structure, while others would be expected to change. Even 
qualitative results can provide important information about the underlying 
physics and strengthen (or challenge) our present theoretical understanding. 
Furthermore, an analysis based on the methods developed here can serve 
as the basis of comparison with numerical codes which simulate the time 
evolution of EP-driven modes. 

Typically, plasma modes cause magnetic, density and temperature fluc¬ 
tuations in the plasma which makes it possible to extensively examine 
them in a well-equipped tokamak such as AUG. Since the observed modes 
cause weak oscillations in the measured quantities, the used diagnostics 
tools were carefully selected by taking into consideration their spatial and 
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temporal resolution and their signal-to-noise ratio. The argument of the 
selection and the description of the applied diagnostic tools are presented in 
section 2.1. The analysis in question requires sophisticated signal processing 
tools. Since this work is focused on the investigation of transient, wave-like 
phenomena, the signal processing methods developed and used are based on 
short time Fourier transform (STFT) [20]. The applied methods are described 
in section 2.2. The amplitude reconstruction method derived in this thesis 
which allowed to investigate the radial structure of the observed modes is 
presented in chapter 3. The dedicated experiments analysed are shown in 
section 4 where the results are also discussed in view of the physical picture. 
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Chapter 2 


Measurement set-up and 
analysis principles 


The frequency of the modes investigated in this thesis is in the order of 
50 - 100 kHz. Both the diagnostic tools and the data processing method were 
chosen to deal with such high frequency oscillations. 

From the diagnostic point of view, fluctuation measurements are required 
which can measure either the magnetic, the density or the temperature 
fluctuations in the plasma with higher temporal resolution than 100 kHz. The 
possible candidates were the magnetic pick-up coil, the reflectrometry, the 
electron cyclotron emission (ECE) and the soft X-ray (SXR) measurements. 
Magnetic measurements do not have spatial resolution, thus it is not suitable 
for the analysis of the radial structure, but it was used to examine the mode 
structure. The reflectrometry diagnostic can measure density fluctuation 
with high temporal resolution, however due to the low spatial resolution 
I only used it to localize the observed modes. The ECE [21] and the 
ECE imaging (ECEI) [22] diagnostics provide local electron temperature 
measurements. These would be good candidates due to their high time- 
and spatial resolution, however BAEs and EGAMs were not visible by ECE 
and ECEI in the investigated shots. Thus, ECE and ECEI are ruled out 
from the analysis. The signal of SXR measurement is a function of electron 
temperature and density. Its sampling frequency is 2 MHz and due to the 
high number of SXR channels it also has a spatial resolution. Therefore, 
this diagnostic was chosen to investigate the changes in the radial structure 
of the modes. The magnetic probes, the soft X-ray measurements and the 
reflectrometry diagnostics are introduced in section 2.1. 

In addition to the suitable diagnostic tools, an appropriate data processing 
method is required. For this purpose, continuous linear time-frequency 
transforms were chosen, because these are ideal to investigate transient 
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signals. The amplitude reconstruction method developed is presented later 
in chapter 3 is based on STFT [20] which is a linear time-frequency 
transform. The mathematical background of these transforms are introduced 
in section 2.2. 


2.1 Measurement setup 

The experiments presented in this thesis were carried out on the ASDEX 
Upgrade tokamak (Axially Symmetric Divertor Experiment Upgrade) [23], 
which is a middle sized tokamak with a major radius of 1.6 m and a 
minor radius of 0.5 m. It started operation at the Max-Planck-Institut 
fiir Plasmaphysik, Garching in 1991 and it plays a substantial role in the 
European fusion program [24], At present, it is Germany’s largest fusion 
experiment. The device is equipped with a flexible neutral beam injection 
(NBI) system which allows to carry out cutting edge research in fast ion 
physics. Furthermore, the tokamak has a wide variety of diagnostics tools. 
Three of them were used in this thesis and these are presented in the 
forthcoming sections. 

2.1.1 Magnetic probes 

The magnetic pick-up coils are mounted on the inner side of the vacuum 
chamber. Depending on the alignment these probes pick up different com¬ 
ponents of the magnetic field fluctuation. On AUG, so-called Mirnov coils 
are placed on the vacuum vessel to measure the poloidal component of the 
magnetic field fluctuations B po \ and so-called ballooning coils are placed closer 
to the plasma on the low field side to measure the radial component B r . Due 
to their relative simplicity and compact size, several coils can be distributed 
over different locations within the plasma vessel, making them suitable for 
the analysis of mode structure. 

The toroidal position of the ballooning coils is illustrated in figure 2.1a 
where a top-down view of the tokamak is shown. These probes are placed in 
almost identical poloidal positions which allows to investigate toroidal mode 
numbers. Furthermore, these coils are placed near the last closed flux surface 
to more efficiently detect the perturbations coming from the plasma. The 
ballooning coils are calibrated by taking into account the transfer functions of 
the different probes [25] . This way the systematic errors in the mode number 
fitting are reduced which allows a better identification of mode numbers with 
reduced number of coils even for transient modes. 
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Figure 2.1: The position of magnetic probes used for the analysis on AUG. 

(a) The six ballooning coils of the toroidal array are indicated with 
red in the top-down view of the tokamak. (b) The 30 Mirnov coils 
indicated with green in a poloidal cross section of the device. Probes 
indicated with empty rectangles were ruled out from the analysis, 
because the investigated modes were not visible in their signal. The 
cyan rectangles indicate the passive stabilizing loops (PSLs). 


In order to determine the poloidal mode number, the poloidal array of 
Mirnov coils was used. Mirnov coils are placed on the vacuum vessel, a little 
further from the plasma and these are not calibrated. In turn, the 30 Mirnov- 
coils of the poloidal array constitute a full poloidal ring around the plasma. 
The position of these probes is shown in figure 2.1b in a poloidal cross section 
of AUG. Probes indicated with empty green rectangles are ruled out from 
the mode number analysis because the observed modes were hardly seen in 
their signals. Most of the ruled out coils have a weak signal because these are 
placed behind the divertor elements or the passive stabilizing loops (PSLs). 
Despite the difficulties concerning the use of Mirnov-coils, due to the high 
number of available probes, the poloidal mode number analysis gives accurate 
results for low mode numbers. 

2.1.2 Soft X-ray measurement 

The magnetic probes are very useful in mode identification and mode number 
analysis, but they cannot provide information on the radial structure of 
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plasma modes. For this purpose, an other fluctuation measurement was 
chosen: the soft X-ray diagnostic system. The major part of the radiation 
coming from a hot fusion plasma is in the soft X-ray range (100 eV to 
10 keV). The energy spectrum of Soft X-ray radiation consists of a continuum 
radiation resulting from Bremsstrahlung, a recombination radiation resulting 
from free-bound transition, and a line radiation resulting from bound-bound 
transition. In a plasma the light intensity from Bremsstrahlung radiation can 
be estimated by the following formula [26]: 

P B r = 1.69 • 10- 32 • ney/feY^ ( Z 2 N(Z )), (2.1) 

where n e is the electron density in 1/cm -3 , T e is the electron temperature 
in eV and Z is the charge number of the given ionization state. The sum is 
executed on all ionization states and the radiation power is given in W/cm -3 . 
Eq. (2.1) tells us that the Bremsstrahlung emitted by a deuterium plasma 
is proportional to the square of the electron density and the square root 
of the electron temperature. However, the soft X-ray radiation also involves 
recombination and line radiation. On AUG, the SXR detectors have a 75 
//m beryllium filter foil in front to suppress low energy photons. This way, 
photon energies below 1 keV are blocked by the beryllium and line radiation 
from light impurities is negligible for the Soft X-ray diagnostic. Thus, they 
affect the Soft X-ray radiation mainly via effective charge number Z c g. On 
the other hand, due to the full tungsten wall coverage of AUG, line radiation 
can even exceed the continuum radiation seen in Soft X-ray diodes. From the 
analysis point of view, it is sufficient to know that the measured radiation is 
a function of electron temperature and density. 

The SXR system of AUG consists of planar photo diode arrays with 
pinhole optics. A diode collects light from a thin cone shaped volume which 
is referred to as the line-of-sight (LOS). AUG is equipped with 8 cameras 
consisting about 200 LOSs [27]. During the analysis of chirping modes I 
examined the phenomenon on the majority of the LOSs that are shown in 
figure 2.2. 

2.1.3 Reflectrometry 

Reflectrometry can give information on the density of the plasma. When an 
electromagnetic wave of a certain frequency propagates through a plasma 
with density increasing in the direction of propagation, it may arrive at a 
point where the electron density equals the cut-off density [8]. At this point 
the wave is reflected and by detecting the reflected wave it is possible to use it 
to diagnose the plasma density and its fluctuation [21], The phase shift of the 
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Soft X-ray Line-of-sights 



Figure 2.2: The line-of-sights (LOSs) of the used soft X-ray channels in the 
poloidal cross-section of AUG. Different colours indicate the different 
cameras. 


reflected wave - which can be measured by using a microwave interferometer 
- gives information on the position of the cut-off layer. 

One operation mode of the diagnostic system is the profile measurement. 
By sweeping the frequency of the incident wave and recording the phase shift 
of the reflected wave as a function of the frequency, a radial density profile 
can be reconstructed. An other operation mode is the density fluctuation 
measurement which is carried out at a fixed frequency. From the time 
variation of the phase shift, density fluctuations near the cut-off layer can 
be determined. If the plasma is probed by several incident waves, the radial 
localization of the mode can be estimated. On AUG there are 4 channels for 
fluctuation measurements both on the low held side and the high held side. 

2.2 Analysis principles 

To handle the rapidly changing mode frequency and the low signal-to-noise 
ratio of the observed EP-driven modes, an advanced amplitude reconstruction 
method was developed. This is based on Short Time Fourier Transform 
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(STFT) which is a linear continuous time-frequency transform. It is ideal 
to investigate the time evolution of transient wave-like phenomena. The 
mathematical background of linear continuous time-frequency transforms 
is presented in section 2.2.1 based on the book of Stephane Mallat: A 
Wavelet Tour of Signal Processing [28]. Then, in section 2.2.2, the mode 
number analysis method which is used to identify EP-driven modes is briefly 
explained. The introduction of analysis principles continues in section 3, 
where the amplitude reconstruction method developed in this thesis is 
presented. 

2.2.1 Short time Fourier transform 

Linear continuous time-frequency transforms are calculated by expanding the 
signal fit) on the basis of families of so-called time-frequency atoms: 

+oo 

Tf{u,i) = (f,g u ,i) = J f(t)g^(t)dt, (2.2) 

— OO 

where g U £ is a time-frequency atom, whose energy is well localized in both 
time and frequency. Variables u and £ are the time and frequency indices 
of the atom identifying its position on the time-frequency plane, and the * 
represents the complex conjugation. The energy density distribution can then 
be calculated by taking the absolute value squared [28]: 

Ef(u,£) = \Tf(u,£)\ 2 . (2.3) 

Ef(u,f) can be interpreted as energy density distribution on the time- 
frequency plane when the signal energy is defined as 

T 

P = J m\ 2 dt (2.4) 

o 

The mode number analysis is based on the phase of the cross-transform which 
is defined in the following way: 

= arg {Tf k (u,f)Tf;(u,f)} , (2.5) 

where fk and fi represent the signals of probes placed in different positions. 

STFT is a type of linear continuous time-frequency transform when the 
family of time-frequency atoms are generated by shifting a real and symmetric 
window g(t) in time (u) and frequency (£): 

9ut{t) = exp {i£t}g(t - u), (2.6) 
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which gives a uniform time-frequency resolution on the time-frequency plane. 
Figure 2.3 illustrates two atoms shifted in frequency (£). Following from 


Re(g(t)e /f 9 Re(g(t)e m ) 



Figure 2.3: Two time-frequency atoms shifted in frequency (£). 


eq. (2.2) STFT can be written as 

+oo 

Sf(u,0 = (/,&*,£> = J f(t)exp{-i£t}g(t 

— OO 


w)df. 


(2.7) 


I applied Gabor atom [28] in the present thesis which means that g(t) is a 
Gaussian function with a a t standard deviation: 


g{t - u) 




exp 


(t ~ u ) 2 ) 
2a? J' 


( 2 . 8 ) 


The energy density distribution defined in eq. (2.3) calculated from STFT 
is called a spectrogram, which is ideal to investigate the time-frequency 
evolution of non-stationary modes [29, 30]. In order to apply a continuous 
transform on sampled time signals, the calculation has to be discretized 
in a way to avoid the degradation of the time-shift invariance property of 
the transform [31]. This is achieved by replacing the integrals with sums 
and discretizing all variables in equation (2.7) with the smallest possible 
steps. The width (2cr t in equation (2.8)) of the Gabor atom determines 
the time-frequency resolution of the transform. According to the specific 
application, the discretization can be done on a more sparse grid, which 
is determined by a parameter called time step, but it has to be at most the 
quarter of the width of the atom to give a continuous transform. Since the 
Gabor atom has an infinite support it has to be truncated [31]. To preserve 
the appealing properties of the continuous transform, only the part where 
its values are several orders of magnitude smaller than the maximum are 
neglected. The resulting length of the window determines the frequency step 
of the transform. In practice the transform is calculated in each time step 
by fast Fourier transform (FFT) using Gaussian window having very long 
support as it is illustrated in figure 2.4. 
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Figure 2.4: The process of calculating STFT on discrete time signals. The 
transform is calculated in each time step by fast Fourier transform 
(FFT) using Gaussian window. 


2.2.2 Mode number analysis 

The evaluation of toroidal and poloidal mode numbers of the observed EP- 
driven modes were carried out by using signals of the magnetic pick-up coil 
system which is described in section 2.1.1. A brief introduction of the mode 
number analysis method 1 used is given in this section, but a more detailed 
description can be found in our recently submitted paper [25]. 

The measurement of magnetic field fluctuations produced by magnetohy¬ 
drodynamic (MHD) modes allows the reconstruction of their mode structure. 
Harmonics of global MHD eigenmodes are generally assumed in the following 
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form [32]: 


Am,n{pfi*,<t>,t) = An,n(p,0*) exp{f(m6>* + 710 ~ Ut)} , (2.9) 

where A is the observable quantity (in our case the magnetic held perturba¬ 
tion) as a function of time t and the so-called straight held line coordinates 
(p,0*,0) [33], A(p,9*) is the radial eigenfuntion and u is the mode frequency. 
The mode structure of the mode is characterized by the m poloidal and n 
toroidal mode numbers. As follows from eq. (2.9), in the case of a single, pure 
sinusoidal, globally coherent mode at hxed ( p,9*,t ), but at different toroidal 
angles, the relative phase between A{4>k) and A (00 is proportional to the 
relative toroidal angle 0 fc — 0 ; of the probe locations and the ratio is the 
toroidal mode number n\ 

vti I = arg {A(0*0} - arg{A(00} = ?i(0 fc - 00 . (2.10) 

I B k= e i 

Similar formula can be constructed for the poloidal mode number at hxed 
toroidal position: 

= arg {A(9 k )} - arg {A(00} = m(9 k - 00 . (2.11) 

Due to its complexity, the 6* transformation of the poloidal coordinates was 
not performed. This causes negligible error, because the investigated modes 
are core localized, thus the effect of the 9* transformation would be small. 

Considering eq. (2.10) and (2.11) the mode number determination can be 
handled as a linear htting problem. In order to evaluate the toroidal mode 
number, one has to determine the phase of the mode in different toroidal 
locations (00, but at the same poloidal position. Eq. (2.10) shows that the 
slope of the linear curve, htted on relative phases between all pairs of signals 
as a function of the relative probe position, gives the toroidal mode number n. 
A similar argument stands for the poloidal mode number by considering 
eq. (2.11). The cross-phases are evaluated from the cross-STFT which is 
defined in eq. (2.5). The method searches for the integer mode number value 
where the residual Q of the fit is minimal: 

Q = J2 \\<Pv~n<fr\\ 2 , (2.12) 

v 

where pp is the phase between pairs of signals, 0-p is the relative probe 
position, n is the toroidal mode number, the sum is executed on the V signal 
pairs and ||.|| is the norm obtained by taking the optimum shift of pp by 
27 tz, where z is an integer number. Eq. (2.12) can also be defined for poloidal 
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mode numbers. The sign of the mode number determines the direction of 
propagation in the device frame, ff the helicity (which is determined by the 
direction of the toroidal magnetic field and the plasma current) is known, the 
sign of the mode number can be related to the ion or electron diamagnetic 
drift direction. By performing the linear fit in each time-frequency point, a 
time-frequency resolved map of the mode numbers can be generated [34], 
Such mode number maps were evaluated to determine the mode numbers of 
the observed EP-driven modes in section 4. 
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Chapter 3 


Amplitude reconstruction of 
chirping waves 


For the analysis of EP-driven modes STFT was chosen because it is ideal to 
examine their time-frequency evolution. In order to investigate the radial 
structure, the mode amplitude has to be determined on different SXR 
LOSs. STFT provides a good basis for the amplitude reconstruction, because 
the time-frequency resolution allows to eliminate the major part of the 
background noise and only the noise component interfering with the mode 
at the mode frequency causes uncertainty in the reconstruction. However, 
to handle the rapidly changing mode frequency and the low signal-to-noise 
ratio, I developed an advanced reconstruction method which is described in 
this section. 

The derivation of the method is explained in section 3.1. A new approach 
for interpolating STFT which is necessary for accurate reconstruction is 
described in section 3.2. In section 3.3, the effect of background noise on 
the reconstruction is discussed. The validation of the method is presented in 
section 3.4, where results on synthetic signals are shown. The linear chirp 
approximation - which is finally applied to investigate chirping EP-driven 
modes - is investigated in more detail in section 3.5. 


3.1 Derivation of the reconstruction method 

In order to perform the reconstruction, a signal model is needed. The 
fluctuation on the soft X-ray signal caused by the observed mode is modelled 
as a frequency and amplitude modulated harmonic wave. The noise of the 
measured signal is dominated by the soft X-ray light coming from the 
background plasma. The noise of the amplifier chain is negligible. The plasma 


background noise is modelled as a Gaussian, additive, white noise z(t). Since 
the application of STFT serves as a narrow band filter, the white noise 
approximation only assumes that the spectral power of the noise is uniform 
in a narrow frequency range (~ 10 kHz). Thus, this is applicable to model 
any kind of broadband noise. The sum of the chirping wave and the additive 
noise takes the following form: 


}(t) = a(t) cos[0(l)] + z(i), 


(3.1) 


where a{t) is the instantaneous amplitude of the wave and the derivative of 
the phase 0(t) gives the instantaneous frequency (d0/dt = cp(t)' = 

During the derivation of the amplitude reconstruction method, the cosine 
function in eq. (3.1) is substituted by exp{i0(t)}. The results can be easily 
converted considering that 

cos = expW(t)}+exp{— t0(t» (3 . 2) 


The reconstruction method is derived in the following way. In eq. (3.1), a(t) 
and 4>(t) are expanded in Taylor series around u: 

“ a W>(u) 

OT ,«(*) = £—D (*-«)'. (3-3) 

U p! 

<hAt) = E - u)>, (3.4) 

VI 

p =0 1 

h,u(t) = a TjU (t) cos[0 T .«(t)], (3.5) 

where (p) denotes the order of the derivation. Then, eq. (3.5) is substituted 
into the definition of STFT (eq. (2.7)). At this point noise z(t) is neglected. 
The effect of the noise on the results will be discussed in section 3.3. By 
evaluating «S7 t )U (w, 0> a W can ex P resse d as a function of the STFT 
transform and the time derivatives of the amplitude and frequency. This 
evaluation was carried out by Oberlin et al. [35, 36] in the first order 
of the frequency. In the present thesis, a more careful treatment of the 
problem is presented. In addition to the expansion of the frequency, the 
time evolution of the amplitude is also taken into account by including its 
Taylor expansion in eq. (3.5). Furthermore, the effect of higher order terms 
is also investigated. Substituting eq. (3.1) into the equation defining STFT 
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(eq. (2.7)) and separating the zeroth order of the amplitude gives: 


Sh,u{u,£) 


+OO 

J a(u ) exp{i</> T>n (f)} exp{— i£t}g(t — u)dt+ 


— OO 




(t ~ u) p exp{i0 T)U (t)} exp{-i£t}g(t 


(3.6) 


u)dt. 


By rearranging this equation, a{u) can be expressed as a function of the 
STFT transform and the time derivatives of the amplitude and frequency: 


a(u) = (3.7) 

7° ~ a (p)( u ) 

Sh, u ( u >£) ~ / , i \ t~u) p exp{i(f> TjU (t)}exp{-i^t}g(t-u)dt 

J p =i P' 

—oo ^ 

+oo 

J exp{i</>T, u (t)} exp{— i£t}g(t — u)dt 

— OO 

Eq. (3.7) gives an exact formula to evaluate the instantaneous amplitude 
if the background noise power is zero. However, this formula requires the 
knowledge of the instantaneous frequency and its derivatives (</> (p> (u)) and the 
derivatives of the amplitude (a^ p \u), p > 1). Therefore, eq. (3.7) is a recursive 
definition of a(u), because the knowledge of the amplitude is required to 
calculate its derivatives. An estimation for the frequency and the amplitude 
can be given considering eq. (3.6) in the zeroth order for the amplitude and 
in the first order for the phase: 

+oo 

Sfr,u{ u ,£) ~ j a(u) exp{i(j)(u)+i(j)'(u)(t — u)} exp{—i£t}g(t — u)dt. (3.8) 

— OO 

This formula gives exact instantaneous amplitude in the case of constant 
amplitude and constant frequency waves. Considering the formula of the g{t) 
Gaussian window function from eq. (2.8) the integral in eq. (3.8) can be 
analytically evaluated: 

SfT,u ( u >0 = a ( u ) \j 2cr tv /7f exp{icp(u) - i£u} exp j -y(0'(tt) - f) 2 j ■ ( 3 - 9 ) 

It is clear from eq. (3.9) that |<S'/t jU (w, £)[ has its maximum where 
Therefore, the frequency of the mode can be estimated by searching for 
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the maxima on the spectrogram. It follows from eq. (3.9) that the mode 
amplitude can be evaluated from the maximum value of |S'/t>(w, £) |: 


, ^ \Sh,u{u,£ = (p'(u))\ 

alu) = - 2 — . -. 

By using this approximation of the amplitude and frequency, the derivatives 
in eq. (3.7) can be estimated and a(u) can be evaluated numerically. 

In principle, the determination of p in eq. (3.7) above which the higher 
order terms can be neglected, depends on the behaviour of the chirping 
wave. The Gaussian window (eq. (2.8)) applied for STFT has a well-localized 
energy in time which time scale is characterized by the oy standard deviation. 
Eq. (3.7) can give a good approximation with a limited p = p max if the 
remaining terms p < p max accurately describe the chirping wave on the time 
scale of the window width a t . For example, if a linear chirp model (first order 
in amplitude and second order in phase) was a good approximation of the 
observed wave on the time scale of oy, then the higher order terms would give 
negligible contribution to the result. 

In summary, the amplitude reconstruction method - if higher order terms 
are taken into consideration - comprises the following steps: First, the instan¬ 
taneous frequency of the mode is estimated by using a maximum-searching 
algorithm. The time evolution of the maxima of the mode frequency is 
traced, which is called the ridge of the STFT transform. Second, by using 
the zeroth order approximation (eq. (3.10)) the instantaneous amplitude is 
estimated from the STFT values of the ridge. Then, the derivatives of the 
amplitude (a^(u)) and frequency (0^(u)) are evaluated from the zeroth 
order approximation. Finally, the instantaneous amplitude is approximated 
by calculating eq. (3.7). In order to avoid the effect of quantization error in 
the numerical derivatives of the amplitude and frequency, an interpolation of 
STFT is required. 



3.2 Interpolation of STFT 

The interpolation method that I used in this thesis could have a broader 
applicability than a simple interpolation algorithm. With this technique the 
reconstruction of the continuous STFT of a signal at any time-frequency 
point is possible. This method is based on a theorem introduced in the book of 
Mallat [28]. Let T be an element of square integrable functions (T € L 2 (M 2 )). 
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There exists / G L 2 (M) such that = Sf(u,£), if and only if, 

OO OO 

^(«o,fo) = ^ I f '&(u,£)K(u 0 ,u,£ 0 ,£)dud€, (3.11) 

— OO —OO 


where 

K(uq,u,£ o,0 = (9u£,9u o,&) ( 3 - 12 ) 

is the reproducing kernel. This theorem shows that not any two dimensional 
distribution can be interpreted as the STFT of some / function. If the 
discretization of the STFT is handled carefully as it is described in 
section 2.2.1, then the integral in eq. (3.11) can be well approximated by the 
discrete STFT. In this way the continuous transform can be reconstructed 
at any given time-frequency point. 

I searched for the local maxima of the absolute value of the STFT at 
each time-point by interpolating the transform in frequency. The method is 
illustrated in figure 3.1. The red circles show the values of the discrete STFT 
and the grey solid line presents the continuous transform. The transform is 
interpolated in the frequency points halfway between the red points. This 
gives the cyan points which are the results of the first iteration of the 
interpolation. Then, the highest value from the cyan points is selected. After 
that the transform is interpolated again halfway between the neighbouring 
cyan points. This gives the green points which determine a more accurate 
estimation for the STFT ridge. In practice, this iteration was repeated 10 — 15 
times to avoid quantization error in the numerical time derivatives of the 
amplitude and frequency. 

3.3 Background noise level 

The derivation presented in section 3.1 is valid if there is no additional noise 
in the signal. Thus, the effect of the noise has to be investigated. As it was 
mentioned before, my signal model (eq. (3.1)) contains a Gaussian, additive, 
white noise: 


f(t) = a(t) cos[0(t)] + z(t) = /chirp (t) + z{t). (3.13) 

Since STFT is a linear transform: 

s (/chirp {t) + z(t)) = 5/chirp(i) + Sz(t). (3.14) 

To illustrate the problem, the transform of the terms in eq. (3.14) are 
calculated separately, then the behaviour of their sum is investigated on 
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Figure 3.1: Illustration of the interpolation method which was used to accurately 
evaluate the STFT ridge. 


the complex plane. At a given time-frequency point (u, £) the value of the 
transform of f r h\r V (t) gives a complex number whose amplitude and phase 
are deterministic. This complex number is plotted as a white vector on the 
complex plane in figure 3.2a. The STFT of the Gaussian, additive, white 
noise z(t) at a given time-frequency point (w,£) is a complex number whose 
amplitude is proportional to the noise power, and its phase is an uniformly 
distributed random number. The distribution of this vector is illustrated 
with the pink spot on the complex plane in figure 3.2a centred around the 
deterministic component. 

Since the STFT of the Gaussian white noise has a 2D Gaussian 
distribution on the complex plane, the magnitude of the sum of the two 
vectors has a Rice distribution [37]. The Rice or Rician distribution is 
the probability distribution of the magnitude of a circular bivariate normal 
random variable with potentially non-zero mean [38]. Its probability density 
function is the following: 


f(x\u,s) = 4 exp 
s z 


-(x 2 +v 2 n, 
2S2 r 



(3.15) 


where u is the distance of the centre of the 2D Gaussian distribution from the 
origin, s is the standard deviation of the 2D Gaussian distribution and Iq{z) is 
the modified Bessel function of the first kind with order zero. The probability 
density of the Rice distribution with parameters v = 5 and s = 3 is shown 
in figure 3.2b. The mean of a random variable x with Rice distribution is 

E[x] = s\fn~/2 Li/ 2 (—v 2 /2s 2 ) (3.16) 
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Figure 3.2: (a) This figure illustrates the sum of the STFT transform of a 
deterministic wave (white vector with length of v = 5) and a 
Gaussian, additive, white noise (pink spot with standard deviation 
g = 3) at a given time-frequency point (tt,£). (b) The probability 
density of the Rice distribution with parameters v = 5 and s = 3. 


and its variance is 

Var[;r] = 2s 2 + v 2 - ^-Ly 2 (-z/ 2 /2s 2 ), (3.17) 

where L q (x) denotes a Laguerre polynomial. 

The background noise level, i.e. the s of the 2D Gaussian distribution is 
determined in the following way. A time-frequency interval is defined on the 
spectrogram of the given signal where no coherent mode appears. Since it is 
assumed that only the Gaussian white noise gives significant contribution to 
the values over the selected time-frequency interval, the real or imaginary part 
of these values has a Gaussian distribution with standard deviation s, where s 
is the standard deviation of the 2D Gaussian distribution of the noise. Thus, 
the background noise level is estimated by calculating the standard deviation 
of the real part of the STFT values over the selected time-frequency interval. 

As is becomes large or s becomes small the mean of the Rice distribution 
becomes is and the variance becomes s 2 [39]. In the investigated real cases, 
the ratio of v and s is around 10 where the oscillation amplitude is maximal. 
Therefore, in the present thesis, I use the is s limit which gives a good 
approximation to estimate the uncertainty of the reconstructed amplitude. 
In figures showing the reconstructed amplitude in section 4, the error bars 
were derived from the standard deviation s of the background noise by taking 
into account the error propagation in the reconstruction formula. 
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3.4 Numerical tests 


I performed numerical tests to investigate the effect of the higher order terms 
and the background noise on the amplitude reconstruction. These tests were 
carried out on synthetic signals which had similar frequency and amplitude 
time evolution as the measured signals. The amplitude and frequency of the 
generated signals presented in this thesis were the following: 

u>(t) = 2 tt ■ 10 4 (exp{ — lCT 3 f} + 7) [Hz], (3.18) 

a(t) = — exp{5 • 10 2 t} + 10. (3.19) 

First, the tests were carried out without adding noise to the signal to 
investigate the effect of higher order terms in the amplitude reconstruction. 
In figure 3.3a the pre-dehned amplitude of the synthetic chirp is shown with 
red solid line. The other lines show the reconstruction taking into account 
increasing number of higher order terms. Maximum order corresponds to 
the highest derivative of the amplitude and frequency taken into account. 
Figure 3.3a shows that the first order approximation gives a better estimation 
of the pre-dehned amplitude evolution. However, the second order terms do 
not give significant contribution to the results. 


Amplitude 


Amplitude 




Time [s] 


Figure 3.3: (a) The result of amplitude reconstruction taking into account higher 
order terms for a synthetic signal without noise, (b) Amplitude 
reconstruction for synthetic signal with additive Gaussian white 
noise. 


In figure 3.3b the effect of additive noise on the reconstruction is 
examined. In this case an additive Gaussian white noise was added to the 
synthetic signal. The signal-to-noise ratio had similar magnitude as in the 
investigated real measurements. For the synthetic signal, this means that 
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the amplitude of the wave was around 8 and the standard deviation of the 
additive Gaussian noise was set to be 15. These tests had already shown that 
the STFT based method effectively filters out the noise and can deal with 
this very poor signal-to-noise ratio. 

The noise causes higher uncertainty in the zeroth order estimation of 
the amplitude and frequency. This uncertainty leads to high errors when 
the derivatives are calculated. As it is visible in figure 3.3b, the first order 
approximation gives a better estimation for the pre-dehned amplitude than 
the zeroth order approximation. However, the effect of the noise is more 
significant in the case of the second order approximation due to the higher 
order derivatives. 

The results of tests on synthetic signals suggest that the first order 
approximation estimates well the amplitude of the synthetic chirp. This is a 
trade-off between the magnitude of the correction given by the higher order 
terms and the uncertainty arising from the derivatives of the noisy zeroth 
order approximation of the amplitude and the frequency. Since the tests have 
shown that the first order approximation is sufficient for my purposes, this 
is investigated in more detail. 


3.5 Linear chirp approximation 

Eq. (3.6) in the first order for the amplitude and second order for the phase 
gives 


Sh,u( u >£) = 


(3.20) 


+oo 

J ( a(u ) + a'{u){t - u)) exp |z^0(w) + - u) + ^ ^ (t - u ) 2 j j 

— OO 


exp {—iff} 


Vv^t 


exp 


(t ~ uf 
2 ° 2 t 


dt. 


Even in case of a linear chirp, the absolute value of the transform has its 
maximum where (j)'{u ) = £ [35]. Evaluating eq. (3.20) and taking into account 
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that <f>'(u) = £ yields 


Sh,n{u,£) = -y=== exp -£u) 


i<tt'{u) 1 Xu _ a2 


(3.21) 


+oo 

a(w) / exp 


— OO 
+OO 


+a'(w) / (t — w) exp 


- 2^2 )(* r dt + 


#"(«) 1 \u x2 


- ^2) (* ~ u ) r dt 


By performing variable substitution in the integral of the second term for 
t — u, it is clearly visible that this integral contains the product of an odd 
and an even function. Thus, the second term equals to zero which shows that 
the a'(u) does not give contribution to the result. The first integral can be 
evaluated analytically, thus a(t) can be expressed as 


1 1 sfJ‘ f ( n I ^ (j 2 

—--exp {~i{(j){u) - £u)}Sf T , u ( u i 0- (3.22) 

20t\J 7T 

Taking into consideration that a(u) is a real number, its value is given by the 
absolute value of eq. (3.22): 



a(u) 


\sh,uM\ 

V / 2oy7^ 



+ 4>"(u) 2 af. 


(3.23) 


With this formula the amplitude of the wave can be reconstructed when 
the linear chirp is a good approximation on the time scale of cr f . This 
formula differs from the zeroth order approximation (eq. (3.10)) in the 
following correction factor: y/l + (j)"(u) 2 crf. This correction depends on the 
derivative of the frequency and the width of the applied Gaussian window 
which determines the time-frequency resolution of STFT. ft is clear that 
the correction factor is lower if the window width is smaller, i.e when 
the time resolution of the transform is better. However, the window width 
cannot be arbitrary small, because the better the time-resolution the worse 
the frequency resolution. Therefore, a trade-off between the time- and 
frequency resolution is needed. In a previous work I used the zeroth order 
approximation, but the time-frequency resolution was set to minimize the 
effect of the neglected correction [40]. In this thesis I used the first order 
approximation, because it allows to apply a wider window function which 
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provides better frequency resolution, thus more accurate reconstruction of 
the STFT ridge. 

Figure 3.4 shows the effect of the time-frequency resolution on the 
reconstruction. To confirm that the linear chirp approximation gives the 
best estimation for the chirping modes investigated in real signals, in this 
case I reconstructed the amplitude of a real oscillation (BAE) with different 
time-frequency resolution. In figure 3.4a, the zeroth order approximation 


Amplitude Amplitude 



Figure 3.4: (a) The result of amplitude reconstruction by using the zeroth 
order approximation (eq. (3.10)) on a real signal, (b) The result 
of amplitude reconstruction by using linear chirp approximation 
(eq. (3.23)) on a real signal. 

(eq. (3.10)) is used. As it is visible, the reconstructed amplitude depends on 
the time-frequency resolution. Figure 3.4b shows the case of reconstruction 
with the linear chirp approximation (eq. (3.23)). In this case the time- 
frequency resolution has a much weaker effect on the reconstruction. 

As it was discussed in this section, with an accurate signal model, the 
amplitude of rapidly chirping waves can be reconstructed. The effect of 
additive noise is handled separately by taking advantage of the linearity of 
the STFT. Numerical tests had shown that the linear chirp approximation 
gives a good approximation for the amplitude of chirping waves which have 
similar parameters as the ones observed in the measurement. Therefore, in 
the following section, where the experimental observations are discussed, the 
first order approximation (eq. (3.23)) is applied. 
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Chapter 4 

Experimental observation of 
EP-driven modes 


The importance of EP-driven modes regarding future burning plasma 
experiments was discussed in section 1. Present section shows a detailed 
analysis to experimentally characterize the behaviour of EP-driven plasma 
modes. The diagnostic tools presented in section 2 and data processing 
methods developed in section 3 allow to examine the rapid changes in the 
radial structure of bursting EP-driven modes during the non-linear chirping 
phase. 

The experiments in question were carried out on ASDEX Upgrade (AUG). 
The scenario in which BAEs and EGAMs were observed was characterized 
by low density (~ 2-10 19 m -3 ) in order to maximize the fast particle pressure 
fraction. The modes were observed in the current ramp-up phase of the 
discharge when off-axis NBI with 93 keV beam energy was applied. The 
ramp-up phase in a typical AUG discharge lasts about 1 sec. A series of 
discharges with strong Alfven activity were investigated. The list of analysed 
discharges with the approximate time and frequency intervals where BAEs 
and EGAMs were observed is presented in table 4.1. The mode type was 
identified mainly by considering the toroidal mode number. EGAMs have 
a toroidal mode number n — 0 and BAEs have a nonzero toroidal mode 
number. 

The main goal of this thesis is to investigate the changes in the radial 
structure of the mode in the non-linear chirping phase. The only fluctuation 
measurement which has spatial resolution and these modes were visible in its 
signal was the SXR diagnostic. Since the SXR diagnostic is a line-integrated 
measurement, it is not straightforward to reconstruct the radial structure 
of the mode, but many line-of-sights (LOSs) are available which makes it 
possible to qualitatively investigate the time evolution of the radial structure. 
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BAEs 

EGAMs 


Time [s] 

Frequency [kHz] 

Time [s] 

Frequency [kHz] 

#25506 

0.45 - 0.60 

50 - 60 

X 

X 

#28881 

0.64 - 0.72 

70 - 80 

0.85 - 1.00 

60 - 80 

#28884 

0.62 - 0.72 

70 - 85 

X 

X 

#28885 

X 

X 

0.75 - 0.95 

70 - 85 

#30383 

X 

X 

0.25 - 0.60 

30 - 60 

#30945 

X 

X 

X 

X 

#30946 

0.45 - 0.60 

40 - 70 

0.60 - 0.70 

40 - 60 

#30950 

0.45 - 0.55 

50 - 60 

0.70 - 0.80 

60 - 80 

#30951 

0.55 - 0.70 

70 - 90 

X 

X 

#30952 

0.40 - 0.45 

60 - 80 

0.50 - 0.70 

50 - 70 

#30953 

X 

X 

0.55 - 0.75 

40 - 60 

#31213 

X 

X 

0.65 - 1.00 

40 - 60 

#31214 

X 

X 

0.45 - 0.80 

60 - 80 

#31215 

X 

X 

0.45 - 0.75 

40 - 60 

#31216 

X 

X 

0.50 - 0.60 

50 - 70 

#31233 

X 

X 

0.80 - 1.00 

60 - 80 

#31234 

X 

X 

0.65 - 0.80 

40 - 50 


Table 4.1: Investigated shots with strong Alfven activity. The approximate time 
and frequency intervals where BAEs and EGAMs were observed are 
presented. 


The identification of the modes was carried out by using the magnetic 
pick-up coils. Magnetic spectrograms are excellent to follow the time 
evolution of the modes and determine mode numbers. Note that in principle 
EGAMs should only be detectable with density fluctuation measurements 
because it is and electrostatic mode, however, due to sideband coupling it is 
clearly visible on magnetic fluctuation measurements as well [9]. 

In general, the signal-to-noise ratio of the SXR measurements was poor. 
In many cases, modes which were clearly visible on the magnetic spectrogram 
were only visible on one or two LOSs of SXR. The strategy was to find cases 
where the mode is observable on at least 3 adjacent LOSs of a particular 
SXR camera. LOSs of SXR cameras F, G, H, I and J (see section 2.1.2) 
were investigated in the discharges listed in table 4.1 by calculating the 
spectrogram of their signal in the time-frequency range where the mode was 
visible on the magnetic spectrogram. In total, about 1000 SXR spectrograms 
were analysed in this process. Finally, 3 cases for BAEs and 5 cases for 
EGAMs were found where the oscillation amplitude on the SXR signals was 
sufficiently high to identify them visually on the SXR spectrograms. 
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4.1 Beta-induced Alfven eigenmodes 

First, the results from the observation of beta-induced Alfven eigenmodes 
(BAEs) are presented. As it is shown in table 4.1, BAEs were observed in 
seven discharges. However, the signal-to-noise ratio of the SXR signals was 
only appropriate for further analysis in shot #28881. In this shot on-axis NBI 
heating was applied early (from t = 0.35 s), then off-axis NBI was injected 
at t — 0.6 s. Shortly after, bursting BAEs appeared at around 80 kHz as it is 
shown on the magnetic spectrogram on figure 4.1. The modes are observed 
to chirp down in frequency by as much as 10 kHz in 2-3 ms. 
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Figure 4.1: Chirping BAEs with decreasing frequency in the range of 70 — 85 
kHz are visible on the magnetic spectrogram. Numbers denote the 
chirps which were investigated in detail. Modes around 90 kHz were 
not identified. 


4.1.1 Mode number analysis of BAEs 

An important step in mode identification is determining the toroidal and 
poloidal mode numbers. The toroidal mode number was determined from 
the signals with correction [25] of the ballooning coils as it was described 
in section 2.2.2. The result of the toroidal mode number evaluation is 
presented in figure 4.2. In figure 4.2a, the time-frequency resolved mode 
number calculation is shown where it is visible that BAEs have n — 1 toroidal 
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mode number. The mode numbers are plotted only in time-frequency points 
where the residual of the fit is lower than the 5 % of the maximum. This 
filter ensures that only well-fitting mode numbers are taken into account. In 
figure 4.2b, the relative phases between all pairs of signals are plotted as a 
function of relative probe position in the time-frequency point marked with 
a red star in figure 4.2a. The best fitting line with n — 1 slope is plotted with 
black solid line. 


Toroidal mode numbers Toroidal mode number fit 

1.0 

0.5 


o.o 

-0.5 

- 1.0 

0.62 0.64 0.66 0.68 o.O 0.2 0.4 0.6 0.8 1.0 

Time [s] Relative probe positions [re] 

Figure 4.2: Toroidal mode number calculation for BAEs. (a) The result of time- 
frequency resolved mode number calculation. The mode numbers are 
plotted only in time-frequency points where the residual of the fit is 
lower than 5 % of the maximum, (b) The relative phases between all 
pairs of signals are plotted as a function of relative probe position in 
the time-frequency point marked with a red star in figure 4.2a. The 
best fitting line with n = 1 slope is plotted with light green solid line. 




In order to determine the poloidal mode number, the poloidal array of 
Mirnov coils was used. Despite the difficulties concerning the use of Mirnov- 
coils (see section 2.1.1), the results presented in figure 4.3 clearly show that 
the poloidal mode number of the investigated BAEs is m = —3. The negative 
sign means that the mode propagates in the ion diamagnetic drift direction. 
Since the Mirnov coils are not calibrated, in figure 4.3a a different filter 
is applied than in the case of toroidal mode numbers: the mode numbers 
are plotted only in time-frequency points where the value of the minimum 
coherence [41] is higher than 0.3 for number of averages of 3. Figure 4.3b 
shows the relative phases between all pairs of signals as a function of relative 
probe position. 
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Poloidal mode number fit 



Figure 4.3: Poloidal mode number calculation for BAEs. (a) The result of time- 
frequency resolved mode number calculation. The mode numbers are 
plotted only in time-frequency points where the value of the minimum 
coherence is higher than 0.3. (b) The relative phases between all pairs 
of signals are plotted as a function of relative probe position in the 
time-frequency point marked with a red star in figure 4.3a. The best 
fitting line with n = —3 slope is plotted with blue solid lines. 


4.1.2 Changes in the radial structure of BAEs 

The radial structure analysis of BAEs was carried out by using the LOSs of 
SXR camera J. These LOSs are shown in figure 4.4b with red lines. Based 
on the signal-to-noise ratio of the SXR signals three chirps were selected 
for further analysis. The soft X-ray spectrogram of channel J54 is shown 
in figure 4.4a. The instantaneous amplitude of the mode on each LOS is 
calculated by using the first order approximation defined in eq. (3.23). Since 
eq. (3.23) contains the derivative of the instantaneous frequency, the precise 
determination of the instantaneous amplitude requires an accurate frequency 
ridge. Therefore, the frequency ridge was evaluated from the magnetic 
spectrogram, because the magnetic signals have a higher signal-to-noise ratio. 
The ridge was traced by a local maximum searching algorithm [40, 42], 
which follows the mode evolution until the mode amplitude falls below the 
background noise level. 

The time evolution of the instantaneous amplitude of chirp ^2 on the 
different SXR LOSs is shown in figure 4.5a. In order the reduce the effect of 
the noise a smoothing is applied by a moving average with boxcar kernel of 0.5 
ms width. The smoothed amplitudes are presented in figure 4.5b. To examine 
the changes in the radial eigenfunctions, I constructed a radial mapping of the 
mode amplitudes. This means that each LOS is labelled with the normalized 
poloidal flux of the magnetic flux surface to which the LOS is tangential. The 
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SXR spectrogram 



Time [s] 


Soft X-ray Line-of-sights 



Figure 4.4: (a) Spectrogram calculated from the signal of SXR LOS J-54. The 
downchirping BAEs are also visible here similarly to the magnetic 
spectrogram (figure 4.1, note the different time axis.), (b) Six LOSs 
of SXR camera J (red lines) where the signal-to-noise ratio was 
sufficiently high for further analysis. 


Reconstructed amplitudes Smoothed amplitudes 




Time [s] Time [s] 


Figure 4.5: The time evolution of the oscillation amplitude of BAE chirp #2 (see 
figure 4.1) on the different SXR LOSs. (a) Reconstructed amplitudes 
by using linear chirp approximation (eq. (3.23)). (b) Amplitudes 
smoothed by a moving average with boxcar kernel of 0.5 ms width. 


mode amplitude is calculated for each LOS and it is plotted as the function of 
the radial coordinate. This way the radial mapping of the mode amplitudes 
can be evaluated at any time instant. Since this work focuses on the relative 
changes in the radial eigenfunctions of the modes, each curve of the radial 
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map was normalized with its integral. This radial map for chirp ^2 is shown 
in figure 4.6b at time instances indicated with pink arrows in figure 4.5. 
Similarly, figure 4.6a and figure 4.6c show the time evolution of the mode 
amplitude distribution for chirps #1 and #3. The error bars in figure 4.5 
and 4.6 are derived from the estimated background noise level of the SXR 
spectrograms as it was described in section 3.3. The results do not show 
significant changes in the radial distribution of the mode amplitudes. This 
suggests that if there is any change in the radial eigenfunction it is smaller 
than the uncertainty of our measurement. This observation is consistent with 
the physical picture that BAE is a normal mode of the plasma and its radial 
structure strongly depends on the background plasma parameters rather than 
on the EP distribution. 
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Figure 4.6: The radial mapping of the oscillation amplitude of BAEs. The results 
do not show significant changes in the radial distribution of the mode 
amplitudes. 
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4.2 Energetic particle driven GAMs 

EGAMs were observed in the same (or similar) discharges (see table 4.1), 
but in different time-frequency intervals as BAEs. I found 5 cases where 
the signal-to-noise ratio was appropriate for further analysis. The strongest 
EGAMs were observed in discharge #31213, where the angle of the NBI 
beam-line was most off-axis (7.13° respect to the horizontal axis). In this 
discharge 3 consecutive EGAMs were found where the signal-to-noise ratio 
of the SXR signals was sufficient. The evaluated fast ion distribution had a 
peak at p po \ ~ 0.45 according to the fast ion D-alpha (FIDA) spectroscopy 
measurements [43, 44], p po \ = T/To is the normalized poloidal flux which 
acts as a flux label in this case. 


4.2.1 Mode number analysis and localization of 
EGAMs 


For mode identification and mode number analysis again the magnetic pick¬ 
up coils were used. The time evolution of the mode frequency was traced using 
the same ridge-following algorithm as for BAEs. The magnetic spectrogram 
from discharge #31213 which shows the investigated 3 consecutive EGAMs 
is presented in figure 4.7. 
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Figure 4.7: Chirping EGAMs with increasing frequency in the range of 45 - 
60 kHz are visible on the magnetic spectrogram. Numbers denote 
the chirps which were investigated in detail. 
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The mode number analysis of EGAMs was carried out in a similar way 
as it was done for the BAEs. The results of the mode number evaluation are 
presented in figure 4.8, where it is visible that all three EGAMs have n — 0 
toroidal and m = —2 poloidal mode number. 
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Figure 4.8: The result of time-frequency resolved mode number calculation. 

(a) The toroidal mode numbers are plotted only in time-frequency 
points where the residual of the fit is lower than the 3 % of the 
maximum, (b) The poloidal mode numbers are plotted only in 
time-frequency points where the value of the minimum coherence 
is higher than 0.3. 


The next issue to be resolved is to determine the radial location of the 
mode. One estimation of the radial location can be given from the SXR 
measurements. Starting from the plasma core, the first SXR LOS on which 
the mode is not visible gives an outer boundary of the mode location. 
Considering the SXR spectrograms the investigated EGAMs are located 
inside the p po \ = 0.4 surface. The radial location of the mode was also 
estimated from the reflectrometry measurements. Four channels measuring 
density fluctuations were used with fixed frequency (24, 33, 36 and 53 GHz) 
from the low field side. The corresponding density values are 0.71, 1.35, 1.61 
and 3.48 • 10 19 1/m 3 . The detected modes are visible on the spectrograms of 
channels with frequencies 33 and 36 GHz. According to the density prohle 
this means that the mode is located around p po i ~ 0.25 — 0.45. Both results 
are consistent with the FIDA measurement from which it is expected that the 
strongest drive occurs closer to the core than where the fast ion distribution 
peaks (ppeak ~ 0.45). 
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4.2.2 Changes in the radial structure of EGAMs 

The time evolution of the radial structure of EGAMs was investigated in 
the same way as it was done for BAEs. Again, the LOSs of SXR camera 
J were chosen. The selected 6 LOSs where the signal-to-noise ratio was 
appropriate are shown in figure 4.9b in a poloidal cross-section of AUG. The 
soft X-ray spectrogram of channel J54 with 3 consecutive EGAMs is shown 
in figure 4.9a. 
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Figure 4.9: (a) Soft X-ray spectrogram with 3 consecutive EGAMs. (b) The 
position of the selected 6 SXR LOSs in a poloidal cross-section of 
AUG. These channels were used for the radial structure analysis of 
EGAMs. 


The amplitude of the mode is evaluated for all LOSs by using the 
formula defined in eq. (3.23). The reconstructed amplitudes of chirp #2 (see 
figure 4.9a) as a function of time are presented in figure 4.10a for channels 
shown in figure 4.9b. In this figure, the amplitudes were smoothed by a 
moving average with boxcar kernel of 1.25 ms width. The uncertainty of the 
result is indicated by the dashed lines. 

It is already visible in figure 4.10 that the time evolution of the amplitude 
is different on the different channels. Since the SXR measurement is not well 
localized, the exact changes in the radial eigenfunction cannot be tracked, 
but its behaviour can be qualitatively described. The amplitudes shown in 
figure 4.10a are normalized to the amplitude at the beginning of the chirp for 
further analysis. The result is shown on figure 4.10b. In order to investigate 
the changes of the radial structure, the radial mapping was evaluated, i.e. 
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Figure 4.10: (a) The time evolution of the oscillation amplitude of chirping 
EGAM #2 (see figure 4.9a) on the different SXR LOSs. The 
amplitudes were smoothed with a moving boxcar kernel of 1.25 ms 
width, (b) The amplitudes shown in figure 4.10a are normalized to 
the amplitude at the beginning of the chirp (indicated with dashed 
cyan line) for further analysis. 


the normalized amplitudes are plotted as a function of the radial position of 
the channels in different time points. 

The LOSs above and under the magnetic axis are handled separately and 
the time evolution of the radial mappings are shown in figure 4.11. On LOSs 
located above the magnetic axis (figure 4.11a-e), a shrinkage of the mode is 
visible, since as time evolves the relative amplitude in the middle channel is 
rising compared to the outer channels. One example for the radial mapping 
calculated from channels under the magnetic axis is presented in figure 4.Ilf. 
Regarding all 5 cases, radial mappings calculated from channels under the 
magnetic axis do not show any significant change. This is most probably due 
to the longer distance from the observation point to the mode of these LOSs 
as it is visible in figure 4.9b. The LOSs located under the magnetic axis 
observe a bigger volume, detecting more light from the background plasma 
and relatively less from the flux surfaces where the mode is. However, the 
shrinkage of the mode is visible in radial mappings calculated from channels 
above the magnetic axis. As it is visible in figure 4.11a-e, this shrinkage was 
significant in all 5 investigated cases. 
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Figure 4.11: The radial mappings of the oscillation amplitude of EGAMs. 

(a-e) 5 cases showing the radial mapping from SXR LOSs located 
above the magnetic axis. Shrinkage of the mode is visible in all cases, 
(f) One case showing the radial mapping from SXR LOSs located 
under the magnetic. Radial mappings calculated from channels 
under the magnetic axis do not show any significant change. 
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4.2.3 Discussion 


In the following section the theoretical explanation of the experimental results 
are presented. The simulation results shown in this part are provided by 
Philipp Lauber [19]. The observed shrinkage in the mode structure of EGAMs 
during the non-linear chirping phase is consistent with our present theoretical 
understanding. The mode coupled to the EGAM which was measured by the 
magnetic pick-up coils has a poloidal mode number m = —2. However, the 
EGAM itself has a poloidal mode number m = 0. The resonance condition 
to drive EGAMs simply follows from eq. (1.8) if n = 0, m = 0 and l = —1 
are substituted: 

u-uj t = 0, (4.1) 

where c o is the mode frequency and ui t is the transit frequency of the 
interacting particles. This EGAM resonance condition shows that the mode 
frequency is equal to the transit frequency of the interacting ions. 

The exact phase space and real space coordinates of the wave-particle 
interaction are determined by several factors. First, the EGAM drive is 
proportional to the velocity phase space gradient in the EP distribution 
function. Second, the strength of the drive is proportional to the EP density. 
In addition, the spatial dependence of the damping is an important factor. 
Since in case of EGAMs, the main damping mechanism is ion Landau 
damping [16], the damping decreases with temperature. Thus, the damping 
increases towards the plasma center while the EP density - in the discharge 
^31213 - had an off-axis peak. As a result, the radial position of the mode is 
expected to be closer to the magnetic axis than the peak in the EP density. 
From the SXR measurements, the radial position of the mode is expected to 
be at p ~ 0.25, thus the EP distribution was investigated at this position. 

The EP distribution at p ~ 0.25 was calculated by the TRANSP 
simulation code [19, 45] at the time instance when EGAMs appeared. The EP 
density (F(A,E)) times its derivative (dF/d A) is plotted in figure 4.12 as a 
function of the pitch angle (A) and energy ( E ). The EP distribution can excite 
the mode when the drive overcomes the damping. This can be reached where 
the drive is highest, i.e. at the peak of the F -dF/d A function which is around 
(A = 0.75, E = 40 keV) this case as it is shown in figure 4.12. The coordinates 
of the peak of the F ■ OF/ <9A function determine the parameters of interacting 
particles which are passing particles. The coloured curves in figure 4.12 show 
the (A, E) coordinates of particles with different transit frequencies (38, 50 
and 61 kHz). The curves corresponding to higher frequencies move from the 
green curve (38 kHz) towards the pink curve (61 kHz) as it is shown by the 
white arrow. 
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Figure 4.12: The contour plot shows the EP density (F(A,E)) times its deriva¬ 
tive {dF/d A) as a function of the pitch angle (A) and energy ( E) 
at p « 0.25 in discharge #31213 at t = 0.83 s, where EGAMs 
were observed. The coloured curves show the (A, E) coordinates of 
particles with different transit frequencies (38, 50 and 61 kHz). (This 
figure was created from data provided by Philipp Lauber [1 9].) 


The mode frequency of the observed EGAMs starts at around 45 kHz and 
increases until 60 kHz as it is visible in figure 4.7. The curve corresponding 
to the initial 45 kHz frequency would intersect the region of the maximum 
peak in figure 4.12. This means that - considering eq. (4.1) - the observed 
mode frequency is consistent with the simulated EP distribution function 
and the velocity phase space coordinates of the interacting particles are 
approximately (A = 0.75, A = 40 keV). The radial coordinate p ~ 0.25 
and the velocity phase space coordinates (A, E ) determine the orbit of the 
interacting particles. The orbit of EPs is illustrated in figure 4.13 where the 
radial coordinate of the particle is plotted as the function of the normalized 
circulation time. The green curve corresponds to the particles having (p = 
0.25, A = 0.60, E = 30 keV) coordinates and it shows that these interacting 
particles have a A p « 0.34 — 0.17 = 0.17 orbit width. The mode structure of 
EGAMs is determined by the parameters of the drive, i.e. the orbit width of 
the interacting particles corresponds to the mode extent. 

The mode frequency increases until 60 kHz as it is visible in figure 4.7. 
This means that as time evolves, the resonance condition changes and 
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Figure 4.13: The orbit of EPs with different velocity phase space coordinates 
(A ,E). The radial coordinate of the particle is plotted as the 
function of the normalized circulation time. (This figure was created 
from data provided by Philipp Lauber [19].) 


shifts to particles with higher transit frequency. Figure 4.12 shows only the 
equilibrium EP distribution, the evolution of the EP density due to the 
non-linear interaction with fast ion cannot be followed. However, a qualitative 
description can be given. The pink curve in figure 4.12 shows the phase space 
coordinates of particles which have 61 kHz transit frequency. During the time 
evolution of the mode the resonance condition moves from (A = 0.75, E = 40 
keV) toward the pink curve. This means that the resonance condition moves 
to the region of passing particles (particles with smaller A). The exact 
time evolution of the resonance condition cannot be determined from the 
equilibrium EP distribution, but most probably the resonance follows the 
ridge in the F -dF/d A function indicated with the white arrow in figure 4.12. 
More passing particles have a narrower orbit width which can explain the 
experimentally observed shrinkage of the radial mode structure. The orbit of 
EPs with coordinates (p = 0.25, A = 0.25, E = 50 keV) is shown in figure 4.13 
with red. It is visible in this figure that the orbit width corresponding to the 
end of the EGAM time evolution (A p & 0.31 — 0.18 = 0.13) is narrower 
with approximately 25 % than the initial orbit width at the beginning of the 
interaction. The consequence of the narrower orbit width is the shrinkage of 
the mode structure, which is consistent with the experimental observations. 
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Chapter 5 
Conclusions 


The understanding of energetic particle (EP) driven plasma modes plays a 
key role regarding future burning plasma experiments. Super-thermal EPs 
in tokamak plasmas can excite various instabilities, and the most important 
transport process of EPs in the plasma core is their interaction with these 
global plasma modes. The non-linear behaviour of the mode amplitude 
and frequency may exhibit a wide range of different behaviours, which 
significantly influences the impact of the instabilities on the fast particle 
transport. Therefore, in order to comprehensively understand the non-linear 
behaviour of EP-driven instabilities, the investigation of these modes is 
essential. In this thesis the rapid changes in the radial structure of beta 
induced Alfven eigenmodes (BAEs) and EP-driven geodesic acoustic modes 
(EGAMs) were experimentally investigated during the non-linear chirping 
phase. 

For the extensive characterization of the modes three diagnostic systems 
were chosen, namely the magnetic pick-up coils, the soft X-ray (SXR) cameras 
and the reflectrometry measurements. In order to deal with the transient 
behaviour of the phenomena, short time Fourier transform was chosen as a 
basis of the data processing tools developed in this thesis. The time evolution 
of the radial structure of EP-driven modes was examined by using the SXR 
measurements. The SXR line-of-sights (LOSs) observe different parts of the 
plasma, which made it possible to analyse the spatial structure of the modes. 

To follow the amplitude of the oscillation caused by the mode on the 
different LOSs, I developed an amplitude reconstruction method which can 
handle the rapidly changing mode frequency. The derivation of the method 
and its validation with synthetic signals was presented in detail. The effect 
of additive white noise on the reconstruction was also discussed in the thesis. 
The uncertainty of the measurements was estimated from the properties of 
the background noise. 
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I investigated BAEs and EGAMs which were observed in the ramp-up 
phase of off-axis NBI heated plasmas in ASDEX Upgrade. In order to identify 
the instabilities, a thorough mode number analysis was carried out. Then the 
time evolution of their radial structure was examined. The radial structure 
analysis showed that in case of the observed downward chirping BAEs the 
changes in the radial eigenfunction were smaller than the uncertainty of the 
measurement. This behaviour is consistent with that the radial structure 
of BAEs — as normal modes - strongly depends on the background plasma 
parameters rather than on the EP distribution. 

In case of rapidly upward chirping EGAMs the analysis consistently shows 
shrinkage of the mode structure. This can be explained by the changing 
resonance condition in the velocity phase space of EPs. The mode structure 
of EGAMs is sensitive to the EP distribution. The rising frequency of the 
mode indicates that, as time evolves, the EGAM is driven by more passing 
particles which have narrower orbit width. This leads to the experimentally 
observed shrinkage of the mode structure. 

The results shown in the thesis will be presented at the 42nd EPS 
Conference on Plasma Physics [46]. My proposal to carry out further 
experiments in order to increase the event statistics of my investigation 
has been accepted. Thus, dedicated discharges will be performed for further 
analysis of the examined chirping EP-driven modes in the subsequent 
campaign of the ASDEX Upgrade tokamak. 
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